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Abstract 

Three types of Fermions play a fundamental role in our understanding of nature: Dirac, Majo- 
rana, and Weyl. While Dirac fermions are known, the latter two have not been observed as any 
fundamental particle in high energy physics and have emerged as a much-sought-out treasure in 
condensed matter physics. A Weyl semimetal is a novel crystal whose low-energy electronic exci¬ 
tations behave as Weyl fermions. It has received worldwide interest and is believed to open the 
next era of condensed matter physics after graphene and three-dimensional topological insulators. 
However, experimental research has been held back because Weyl semimetals are extremely rare 
in nature. Here, we present the experimental discovery of the Weyl semimetal state in niobium ar¬ 
senide (NbAs). Utilizing the combination of soft X-ray and ultraviolet photoemission spectroscopy, 
we systematically study both the surface and bulk electronic structure of NbAs. We experimen¬ 
tally observe both the Weyl cones and Weyl nodes in the bulk and the Fermi arcs on the surface of 
this system. Our ARPES data, in agreement with our previous theoretical prediction and present 
band structure calculations, provide conclusive evidence for the topological Weyl semimetal state 
in NbAs. Our discovery not only paves the way for the many fascinating topological quantum phe¬ 
nomena predicted in Weyl semimetals, but also establishes a new cornerstone in the correspondence 
between high-energy and condensed matter physics. 
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Weyl semimetals have received significant attention in recent years because they extend 
the classification of symmetry-protected topological phases beyond insulators, demonstrate 
unusual transport phenomena and provide an emergent condensed matter realization of 
Weyl fermions, which do not exist as fundamental particles in the standard model [l? - 
211]. Such kind of topologically nontrivial semimetals are believed to open a new era in 


condensed matter physics. In contrast to topological insulators where only the surface 
states are interesting, a Weyl semimetal features unusual band structure in the bulk and 
on the surface, leading to novel phenomena and potential applications. This opens up 
unparalleled research opportunities, where both bulk and surface sensitive experimental 
probes can measure the topological nature and detect quantum phenomena. In the bulk, a 
Weyl semimetal has a band structure with band crossings, Weyl nodes, which are associated 
with definite chiral charges. Unlike the two-dimensional Dirac points in graphene, the surface 
state Dirac point of a three-dimensional topological insulator or the three-dimensional Dirac 
points in the bulk of a Dirac semimetal, the degeneracy associated with a Weyl node does not 
require any symmetry for its protection other than the translation symmetry of the crystal 
lattice. The low-energy quasiparticle excitations of a Weyl semimetal are chiral fermions 
described by the Weyl equation, well-known in high energy physics, which gives rise to a 
condensed matter analog of the chiral anomaly associated with a negative magnetoresistance 


in transport 


[l7? 


2l| . On the surface, the nontrivial topology guarantees the existence of 


surface states in the form of “Fermi arcs”, which are open curves that connect the bulk 
Weyl nodes. The Fermi arcs on the surface of a Weyl semimetal are by themselves of great 
interest, because Fermi surfaces have to be closed contours in any purely two-dimensional 
band structure. Moreover, they are further interesting because they lead to novel spin 


po 


arization textures 


22| . unusual quantum interference effects in tunneling spectroscopies 


24l |. and a new type of quantum oscillation 23(|, where electrons move in real space between 


different surfaces of a bulk sample when executing a constant-energy orbit in momentum 
space under an external magnetic field. 

Despite the tremendous interest in Weyl semimetals, experimental research devoted to 
its discovery has been held back due to their rareness in nature. The transition metal 


monoarsenides have recently been suggested as Weyl semimetals 


28-|m|. 


By combining 


low-energy (vacuum ultraviolet) and soft X-ray angle-resolved photoemission spectroscopies 
(ARPES), we directly observe both the Weyl cones in the bulk and the Fermi arc surface 
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states on the (001) cleaving surface of NbAs. Specifically, our soft X-ray data clearly reveals 
that the bulk bands touch in discrete points in the bulk Brillouin zone (BZ) and disperse 
linearly along both the in-plane and out-of-plane directions, which, therefore, form the Weyl 
cones. Our low-energy ARPES measurements reveal the existence of Fermi arc surface states 
and show that they connect the projected bulk Weyl nodes. Our systematic ARPES data, 
and its excellent agreement with our previous theoretical prediction and hrst-principle band 
structure calculations, collectively demonstrate the Weyl semimetal state in NbAs. Our 
discovery not only paves the way for the many fascinating topological quantum phenom¬ 
ena predicted in Weyl semimetal state 17H22|. but also establishes a new keystone in the 
correspondence and synergy between high-energy and condensed matter physics. 

Niobium arsenide, NbAs, crystallizes in a body-centered tetragonal Bravais lattice, space 
group lAimd (109), point group € 4 ^. Our X-ray diffraction (XRD) obtains lattice constants 
of a = 3.45 A and c = 11.68 A, consistent with two earlier crystallographic studies 


26 


271]. The basis consists of two Nb atoms and two As atoms (Fig. [T^). The lattice can be 


understood as a stack of square lattice layers, where each layer consists entirely of either 
Nb or As atoms. In the crystal, each layer is shifted by half a lattice constant in either 
the X OT y direction relative to the layer below it. These shifts give the lattice a screw 
pattern along the z direction, which leads to a non-symmorphic C 4 rotation symmetry that 
requires a translation along the 2 ; direction by c/4, or one-fourth a lattice constant in the 
conventional unit cell. More importantly, it can be seen that the NbAs crystal lacks space 
inversion symmetry. Since the system respects time-reversal symmetry, the broken inversion 
symmetry is a crucial condition, without which the Weyl semimetal phase is not possible. 
Fig. [TJa shows a scanning tunneling microscope image on the cleaved (001) surface. The 
topography image clearly resolves a square lattice without any obvious defect, demonstrating 
the high quality of our sample. From the topography, we obtain a lattice constant a = 3.52 
A, consistent with the value determined by XRD. 

We present an overview of the band structure of NbAs before discussing our ARPES 
results. Fig. [T]c shows the calculated band structure without spin-orbit coupling. It can be 
seen that the conduction and valence bands cross each other near the E, S' and N points, 
in agreement with the semimetal groundstate. In order to understand the detailed conhg- 
uration of the band crossing in momentum space, we have calculated the band structure 
throughout the bulk BZ. As shown in the left panel of Fig. [T]d, without spin-orbit coupling. 
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the conduction and valence bands dip into each other, giving rise to four ID ring-like cross¬ 
ings at the kx = Q and ky = Q planes. After the inclusion of spin-orbit coupling (The left 
panel of Fig. [Ud), the bands become gapped everywhere along the rings, but a number of 
point-like band crossings, the Weyl nodes, emerge in the vicinity of the rings. Our calcula¬ 
tion shows that there are in total 24 Weyl nodes, which can be categorized into two groups. 
We name the 8 Weyl nodes that are located on the kz = ‘I'njc plane as Wl, and the other 
16 away from the kz = ‘It: jc plane as W2. We further study the projection of these Weyl 
nodes onto the (001) surface. All Wl Weyl nodes project as single Weyl nodes in the close 
vicinity of the surface BZ edge X and Y points. On the other hand, interestingly, each pair 
of W2 Weyl nodes with the same chiral charge are projected onto the same location on the 
(001) surface BZ. Therefore, there are 8 projected W2 Weyl nodes with the projected chiral 
charge of ±2. The distribution of the projected Weyl nodes and their chiral charges are 
shown in the right panel of Fig. [T)d. The theoretically calculated and ARPES (001) Fermi 
surface measured at the incident photon energy of 50 eV are shown in Figs. [1^, respectively. 
A good agreement is found between the ARPES and the calculation results. We discuss 
some general features of the Fermi surface before going into details. The surface BZ is found 
to be a square. Moreover, the momentum space distance between the surface BZ center F 
and the edge X or F is about 0.91 A ^ in our ARPES data. These observations are in 
agreement with the tetragonal lattice {a = b = 3.45 A) of NbAs and the lattice constant 
(vr/a = 7r/3.45 = 0.91). Both our ARPES data and calculation show that the Fermi surfaces 
are located near the surface BZ edge X and Y points and the midpoints between the F — X 
and the F — F lines. Furthermore, it can be seen from our data and calculation that the 
Fermi surface violates the C 4 rotational symmetry, even though the surface BZ is a square. 
This is entirely consistent with the crystal structure, where the rotational symmetry is im¬ 
plemented as a screw axis that sends the crystal back into itself after a C 4 rotation and a 
translation by (0, |, |) along the rotation axis. As a result, the (001) surface breaks the 
rotational symmetry of the crystal and the surface Fermi surface does not need to respect 
the C 4 symmetry. 

We now present our ARPES data to demonstrate the existence of the bulk Weyl cones 
and the surface Fermi arcs in NbAs. We have systematically studied the electronic structure 
of NbAs using both low-energy (vacuum ultraviolet, 35 eV to 90 eV in our experiments) and 
soft X-ray (350 eV to 1000 eV in our experiments) ARPES. We found that at the vacuum 
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ultraviolet incident photon energies (e.g. Fig. [T^), the surface states dominate the ARPES 
spectral weight, whereas the bulk bands are very weak. On the other hand, at soft X-ray 
energies the situation is reversed and, therefore, the bulk bands become predominant, which 
is consistent with the enhanced bulk sensitivity at the soft X-ray energies. 

We start by presenting the measurement of the bulk Weyl cones using soft X-ray ARPES. 
Fig. [2^ shows the calculated Fermi surface contour of the bulk bands at the ky = 0 

plane. The Fermi surface consists of identical contours centered the S and S points. We 
note that these contours arise from the ID nodal line crossings at the ky = 0 plane without 
spin-orbit coupling (Fig. [Dd), which explains why their shape is very similar to the nodal 
line crossings. As discussed above in Fig. [T]d, we do not expect to observe the Weyl nodes 
in Fig. [2h because all Weyl nodes are away from the ky = 0 plane. Specifically, W1 nodes 
are located at the ky = 0.005 7 i/a plane whereas W2 nodes are at the ky = 0.012 n/a plane 
according to our calculation. Fig. [JJa shows the ARPES measured kz — k^ Fermi surface 
over multiple BZs along the kz direction. We find a good agreement between our ARPES 
data and the calculation. Furthermore, we emphasize that the clear dispersion along the kz 
direction observed in Fig. |2]b hrmly shows that our soft X-ray ARPES primarily measures 
the bulk bands instead of the surface states. This is accurate because if we were measuring 
the surface states, then one would expect the bands to be kz independent, which means that 
the Fermi surface consists of straight lines that are parallel to the kz axis. 

In our calculation, the location of the W2 is {kx,ky,kz) = (0.56l7r/a, 0.0127r/a, I.IGtt/c). 
We use that as a guideline, and choose an incident photon energy (i.e. a kz value) that 
corresponds to the W2 Weyl nodes. Fig. [2]d shows the ARPES measured kx — ky Fermi 
surface at the kz value that corresponds to the W2 Weyl nodes using incident photon energy 
of 651 eV. We compare the soft X-ray and low energy Fermi surface Fig. |2]d and Fig. [1^. 
Although they are measured on the same (001) surface of the same sample, the soft X-ray 
ARPES Fermi surface contains only Fermi points and is (^4 symmetric and, whereas, the low 
energy ARPES Fermi surface contains many pockets and violates the C 4 symmetry. This, 
again, provides compelling evidence for the soft X-ray ARPES measuring the bulk bands 
while low energy ARPES is sensitive to the surface states. We show the schematic of the 
Fermi surface at the kz that corresponds to W2 nodes in Fig. |2^. Apart from the eight W2 
Weyl nodes, there are additional trivial pockets along the kx and ky axes, whcih are denoted 
by the red circles. This is consistent with our ARPES data in Fig. |2]d. We therefore study 
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the energy dispersion of the W2 Weyl cones by choosing a momentnm space cnt direction, 
Cut 1, which goes across a pair of W2 Weyl nodes, as noted in Fig. [2]d. The dispersion map 
is shown in Fig. [2]f over a wide energy range, where a linear dispersion is clearly observed. 
Since Cut 1 goes through two nearby W2 Weyl nodes, one would expect to observe two 
linearly dispersing cones in the ARPES spectra. However, the two cones are very close to 
each other in momentum space. According to our calculation, the k space distance between 
the two W2 nodes is 0.024 tt/o ~ 0.022 A which is challenging to resolve with soft X- 
ray ARPES. However, due to the sharp ARPES spectrum collected from our high quality 
samples, we are able to barely see that the ARPES dispersion in Fig. |2^ are consistent with 
the scenario with two cones. We further perform a second derivative image of the dispersion 
in Fig. [2]h, where the two cones are better seen. The ky separation of the two Weyl nodes 
is about 0.025 A which agrees with the calculation. Since a Weyl cone features linear 
dispersion along all three directions, it is therefore also important to study the dispersion of 
the W2 Weyl bands along the out-of-plane kz direction. We £x the kx and ky at the location 
that corresponds to a W2 Weyl node as shown by Cut 2 in Fig. [2)d and study the energy 
dispersion as a function of kz- As shown in Fig. 13, the Weyl band is found to be strongly 
kz dispersive, which demonstrates its 3D nature. Moreover, throughout the kz range shown 
in Fig. 13, the band only crosses the Fermi level at two discrete kz, where the kz values of 
these two points are approximately (52 ± 1.2)7r/c = ±1.27r/c. This is again consistent with 
the kz value of the W2 Weyl nodes from calculation, kz = ±1.167r/c. In the vicinity of the 
two kz values where the bands cross the Fermi level, the bands disperse linearly within our 
instrumental kz resolution (Fig. [3)- Our ARPES measured E — kz dispersion data in Fig. 3 
is well consistent with the corresponding theoretical calculation in Fig. [3- 

Next, we use the calculated momentum space location of Weyl nodes Wl, {kx, ky, kz) = 
(0.9767r/a, O.OOdvr/a, Ovr/c) as a guideline and choose a photon energy that corresponds to 
Wl. Fig. |3t shows the kx — ky Fermi surface at the kz value that corresponds to Wl. The 
schematics of the expected Fermi surface at kz =W1 is shown in Fig. 3- It can be seen 
that the Fermi surface at kz =W1 consists of both the Wl Weyl nodes that are very close 
to the surface BZ edge [{kx, ky) = (tt/o, 0) or (0, vr/a)] and some trivial pockets closer to the 
BZ center. Therefore, we study the energy dispersion of the Wl Weyl cone by choosing a 
momentum space cut direction. Cut 3, as noted in Fig. 3k. We observe a clear dispersion 
of the Wl Weyl cone in Figs. 3^1,11. Since Cut 3 goes across two nearby Wl Weyl nodes. 


One should, in principle, observe a pair of Weyl cones in the dispersion map depicted by 
Figs. llJm,!!. However, because the two W1 Weyl nodes are extremely close to each other 
(the k space separation is 0.007 A ^ in our calculation, which is much smaller than in W 2 ), 
our experimental resolution is not sufficient. Therefore, the separation of the two cones for 
W1 cannot be resolved. Moreover, our data shows that the Fermi level is about 38 meV 
above the energy of the W1 Weyl node and, as a result, a small part of the upper Weyl cone 
is seen. This is in agreement with the results from our hrst principle calculation, shown in 
Fig. [2]c;. More precisely, our calculation shows that the energy difference between the W1 
and W2 Weyl nodes is found to be 36 meV. Consequently, the systematic measurements 
presented here, when corroborated by our calculations, experimentally reveal the bulk Weyl 
cones in NbAs. 

After establishing the existence of the bulk Weyl cones, we study the Fermi arc surface 
states in NbAs by using low energy (vacuum ultraviolet) ARPES in Fig. [3l Fig. [3^ shows 
the constant energy contours of the ( 001 ) surface at different binding energies with the 
incident photon energy hxed at 50 eV. It can be seen that the constant energy contours, 
especially the pockets centered at the BZ edge X and Y points, violate the C 4 symmetry, 
which serves as the signature for these contours arising from the surface states. At the 
Fermi level {E-q = 0 meV), our data shows that the Fermi surface consists of three types 
of pockets: a dog-bone shaped contour centered at the Y point, a long elliptical contour 
centered at the X point, and a tadpole-shaped pocket that goes along each the F — X or 
F — y) line from the midpoint of the line toward the X (R) point. Our ARPES data and 
calculation (Fig. H^) show agreement on the tadpole and the dog-bone surface states, but 
not on the elliptical contour. Note, the calculation result can be tuned by adjusting the 
surface potential and other surface parameters to further improve its consistency with the 
data. We focus on the tadpole surface states where the agreement between the data and 
the calculations is the best. As shown in Fig. [3]e, the tadpole surface states give rise to the 
Fermi arcs that are connected to the W2 Weyl nodes. Specihcally, the pair of W2 Weyl 
nodes are located at the joint between the “head” and the “tail” of the tadpole. Therefore, 
each W2 is connected to two Fermi arcs, consistent with its projected chiral charge of ±2. 
Because we have resolved both the W2 Weyl nodes using soft X-ray ARPES (Fig. |2]d) and 
the tadpole surface states using low energy ARPES (Fig. |3h), we, therefore, superimpose 
these two Fermi surfaces together to scale, as shown in Fig. [3b. Indeed, we hnd that the 
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W2 nodes are located between the “head” and the ’’tail” of the tadpole. This demonstrates 
that the observed tadpole surface states consists of three Fermi arcs, namely a curve that 
directly connects the pair of W2 nodes, and two nearly straight lines that runs from each 
W2 toward the surface BZ edge. On the other hand, resolving the Fermi arc connectivity 
associated with the W1 Weyl nodes is much more challenging because the W1 nodes are too 
close to each other (the k space separation is 0.007 A ^). In fact, it is even difficult to figure 
out the connectivity associated with the W1 in calculation, because the W1 nodes are too 
close to each other and because there are additional bulk bands at the /c-space between the 
pair of W1 nodes. Nevertheless, our observation of the tadpole Fermi arc surface states that 
connect the W2 Weyl nodes is sufficient to experimentally establish the existence of Fermi 
arc surface states in NbAs. 

Our direct experimental observation by ARPES of the unusual way in which Fermi arcs 
connect Weyl points on the (001) surface of NbAs places strong constraints on the dispersion 
of the Fermi arcs. To understand how, we recall hrst that the Weyl points W2 are arranged 
in such a way that, for an NbAs sample cleaved on the (001) surface, two Weyl points of 
the same chirality project onto the same point of the surface BZ. As a result, each of the 
Weyl point projections has a chiral charge of ±2. This chirality requires not only that two 
Fermi arcs connect to each of the Weyl point projections in the surface BZ, but also that the 
Fermi arcs connecting to a given Weyl point projection be co-propagating. To clarify this 
requirement, consider a circular pipe passing through the bulk BZ of NbAs which encloses 
two Weyl points of chiral charge +1, as shown in the bottom panel of Fig. 4a. The Chern 
number on this manifold is +2. By cleaving the NbAs sample on the (001) surface, we 
introduce a boundary on this manifold corresponding to a circle in the surface BZ, as shown 
in the middle panel of Fig. 4a. On this circle, the Chern number of +2 requires that the band 
structure have two gapless co-propagating chiral boundary modes, illustrated by the red and 
blue curves in the top panel of Fig. 4a. By determining the chirality of the Fermi arc surface 
states using this circle in the surface BZ, we see that if the Weyl points are connected by 
tadpole Fermi arcs, as we observe in our ARPES spectra, then the Fermi arcs must disperse 
as illustrated in Fig. 4b. The green cones correspond to the Weyl cones projected onto the 
surface BZ, while the red and blue sheets correspond to the tadpole Fermi arcs. We can 
better visualize the evolution as a function of binding energy by considering cuts through 
this band structure. We illustrate in Fig. 4c constant-energy contours of this band structure 
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at a cut above the energy of the Weyl points e > ew, near the energy of the Weyl points 
e ~ Ew and below the energy of the Weyl points e < Ew We hnd that the condition for 
co-propagating chiral modes actually requires the red surface states to grow with binding 
energy, while the blue surface states shrink with binding energy. Our prediction of this 
structure of the tadpole Fermi arc surface states is consistent with our numerical calculation 
of the chiral charge of Weyl points in NbAs and our direct experimental observation by 
ARPES of the unusual way in which Fermi arcs connect Weyl points on the (001) surface 
of NbAs. 

ACKNOWLEDGEMENT 

MZH, SYX and IB acknowledge theoretical discussions with I. Klebanov, A. Polyakov, 
P. Steinhardt and H. Verlinde. We gratefully thank Sung-kwan Mo, Jonathan Denlinger, 
for beamline support at the ALS beamlines 10.0.1 and 4.0.3, respectively. We acknowledge 
C.-H. Hsu for technical assistance with some parts of the theoretical calculations presented 
here. Grants: 


[1] Weyl, H. Elektron und gravitation. I. Z. Phys. 56, 330-352 (1929). 

[2] Nielsen, H. B. &: Ninomiya, M. The Adler-Bell-Jackiw anomaly and Weyl fermions in a crystal. 
Phys. Lett. B 130, 389 (1983). 

[3] Balents, L. Weyl electrons kiss. Physics 4, 36 (2011). 

[4] Wilczek, F. Why are there analogies between condensed matter and particle theory? Phys. 
Today 51, 11 (1998). 

[5] Turner, A. M. & Vishwanath, A. Beyond band insulators; topology of semi-metals and inter¬ 
acting phases. arXiv:1301.0330 (2013). 

[6] Haldane, F. D. M. Attachment of surface “Fermi arcs” to the bulk Fermi surface; “Fermi-level 
plumbing” in topological metals. arXiv:1401.0529 (2014). 

[7] Xu, S.-Y. et al. Observation of “Fermi arcs” surface states in a topological metal. Science 347, 
294-298 (2015). 

[8] Hasan, M. Z., Xu, S.-Y., Neupane, M. Topological insulators, topological crystalline insulators, 



11 


topological Kondo insulators, and topological semimetals. arXiv:1406.1040 (2014). 

[9] Wan, X., Turner, A. M., Vishwanath, A., Savrasov, S. Y. Topological semimetal and Fermi- 
arc surface states in the electronic structure of pyrochlore iridates. Phys. Rev. B 83, 205101 
( 2011 ). 

[10] Murakami, S. Phase transition between the quantum spin Hall and insulator phases in 3D: 
emergence of a topological gapless phase. New J. Phys. 9, 356 (2007). 

[11] Burkov, A. A. &: Balents, L. Weyl semimetal in a topological insulator multilayer. Phys. Rev. 
Lett. 107, 127205 (2011). 

[12] Xu, G. et al. Chern semi-metal and quantized anomalous Hall effect in HgCr 2 Se 4 . Phys. Rev. 
Lett. 107, 186806 (2011). 

[13] Singh, B., Sharma, A., Lin, H., Hasan, M. Z., Prasad, R., &: Bansil, A. Topological electronic 
structure and Weyl semimetal in the TlBiSe 2 class of semiconductors. Phys. Rev. B 86, 115208 
( 2012 ). 

[14] Bulmash, D., Liu, C.-X., & Qi, X.-L. Prediction of a Weyl semimetal in Hgi-j^-yCdjjMn^Te. 
Phys. Rev. B 89, 081106(R) (2014). 

[15] Liu, J. &: Vanderbilt, D. Weyl semimetals from noncentrosymmetric topological insulators. 
Phys. Rev. B 90, 155316 (2014). 

[16] Xu, S.-Y. et al. Topological phase transition and texture inversion in a tunable topological 
insulator, TlBiSe 2 . Scienee 332, 560-564 (2011). 

[17] Zyuzin, A. A. Sz Burkov, A. A. Topological response in Weyl semimetals and the chiral 
anomaly. Phys. Rev. B 86, 115133 (2012). 

[18] Wei, H. , Chao, S.-P., &: Aji, V. Excitonic phases from Weyl semimetals. Phys. Rev. Lett. 109, 
196403 (2012). 

[19] Parameswaran, S. A. et al. Probing the chiral anomaly with nonlocal transport in three- 
dimensional topological semimetals. Phys. Rev. X 4, 031035 (2014). 

[20] Ashby, P. E. C. &: Carbotte, J. P. Magneto-optical conductivity of Weyl semimetals. Phys. 
Rev. B 87, 245131 (2013). 

[21] Vazifeh, M. M. &: Franz, M. Electromagnetic response of Weyl semimetals. Phys. Rev. Lett. 
Ill, 027201 (2013). 

[22] Ojanen, T. Helical Eermi arcs and surface states in time-reversal invariant Weyl semimetals. 
Phys. Rev. B 87, 245112 (2013). 



12 


[23] Potter, A. C., Kimchi, L, & Vishwanath, A. Quantum oscillations from surface Fermi-arcs in 
Weyl and Dirac semi-metals. arXiv: 1402.6342 (2014). 

[24] Hosur, P., Friedel oscillations due to Fermi arcs in Weyl semimetals. Phys. Rev. B 86, 195102 

( 2012 ). 

[25] Cho, G. Y., Bardarson, J. H., Lu, Y.-M., & Moore, J. E. Superconductivity of doped Weyl 
semimetals; finite momentum pairing and electronic analog of the ^He-A phase. Phys. Rev. B 
86, 214514 (2012). 

[26] Furuseth, S. & Kjekshuh, A. The crystal structure of NbAs. Acta Crystallogr. 17, 1077-1078 
(1964). 

[27] Boiler, H. & Parthe, E. The transposition structure of NbAs and of similar monophosphides 
and arsenides of niobium and tantalum. Acta Crystallogr. 16, 1095-1101 (1963). 

[28] Huang, S.-M. et al. An inversion breaking Weyl semimetal state in the TaAs material class. 
arXiv: 1501.00755 (2015). 

[29] Weng, H. et al. Weyl semimetal phase in non-centrosymmetric transition metal monophos¬ 
phides. arXiv:1501.00060 (2015). 

[30] Xu, S. Y. et al. Experimental realization (Discovery) of a topological Weyl semimetal phase 
with Eermi arc surface states in TaAs. arXiv: 1502.03807 (2015). 

[31] Lv, B. Q. et al. Discovery of Weyl semimetal TaAs. arXiv:1502.04684 (2015). 

[32] Xu, S.-Y. et al. Observation (Discovery) of Eermi arc surface states in a topological metal. 
Science 347, 294-298 (2015). 

[33] Hasan, M. Z. & Kane, C. L. Topological Insulators. Rev. Mod. Phys. 82, 3045-3067 (2010). 



13 


a 14^md (ooi) view 5 STM topography 




e NbAs Fermi arcs 



FIG. 1: Topological electronic structure of NbAs: Weyl nodes and Fermi arcs, a, 

Body-centred tetragonal structure of NbAs, shown as stacks of Nb and As layers, b, Scanning 
tunneling microscopy (STM) topographic images of cleaved surfaces, c, First-principles band 
structure calculation of the bulk NbAs without spin-orbit coupling, d, (Left) schematics of the 
distribution of the Weyl nodes in the three-dimensional Brillouin zone (BZ) of NbAs. The red and 
blue lines represent the nodal lines without considering spin-orbit coupling. (Right) schematics 
showing the projected Weyl nodes and their chiral charges on the 
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FIG. 1: (001) Fermi surface of NbAs. The projected Weyl nodes are denoted by black and white 
circles whose color indicates the opposite chiral charges of the Weyl nodes, e, First-principles 
band structure calculated (left) and the ARPES measured Fermi surface (right) of the (001) Fermi 
surface of NbAs. The Fermi arcs are clearly resolved in our ARPES measurements in agreement 
with the theoretical prediction. 
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FIG. 2: Weyl cones in NbAs. a, First-principles calculated Fermi surface map at = 0 in the 
kz-kx plane, b, Soft X-ray ARPES Fermi surface map at /cy = 0 in the kz-kx plane, which agrees 
well with the theoretical prediction in a. The measurements were conducted 
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FIG. 2: using photon energies (oc kz) of 518 - 800 eV. c, First-principles calculations of the energy- 
momentum dispersions of the two types of Weyl nodes (W1 and W2) in NbAs which show that 
these nodes are offset in energy by 36 meV relative to each other, d, Soft X-ray ARPES Fermi 
surface map at kz = W2 in the kx-ky plane, revealing the locations of the W2 Weyl nodes, e, 
Schematics corresponding to d indicating the locations of the W2 Weyl nodes (black and white 
circles marked as “2”) and other trivial bulk bands (red dots), f, Soft X-ray ARPES spectra, 
g, its zoomed-in version close to the Eermi level and h, its curvature plot along Cut 1 direction 
shown in d that goes through the twin Weyl nodes marked as “2” in e. The photon energy of the 
measurements in d, f, and g is 651 eV. i, Soft X-ray ARPES spectra along Cut 2 direction shown in 
d (normal to the plane), which again reveals the existence of two Weyl cones along the kz direction, 
j, Eirst-principles calculations of the energy-momentum dispersions of the W2 Weyl nodes along 
kz corresponding to the ARPES measurements in i. k, Soft X-ray ARPES Eermi surface map at 
kz = W1 in the kx-ky plane, revealing the locations of the W1 Weyl nodes. The photon energy 
here is 611 eV. 1, Schematics corresponding to k indicating the locations of the W1 Weyl nodes 
(black and white circles marked as “1”) and other trivial bulk bands (red dots), m, Soft X-ray 
ARPES spectra and n, its zoomed-in version close to the Eermi level along Cut 3 direction shown 
in k that goes through the Weyl nodes marked as “1” in 1. As can be seen in n the energy of the 
W1 Weyl node is below the Eermi level and W2, in agreement with the calculations in c. 



17 



b Weyl nodes and Fermi arcs c Unusua^^c Topology 




^ 0.5 


0 



0.4 




0.8 -I 


-1 -0.5 1) oY 

Tadpole Fermi arcs along T - X 
e _f 


^Arc 

• •^W2 



1^1 


W2 


X 


- 0.2 0 „ 0.2 



^-space monopoles 
& antimonopoles 



-0.5 0.0 o 0.5 1.0 

T(A-Y) 


FIG. 3: Observation of Fermi arc surface states on the (001) surface of NbAs. a, High- 
resolution ARPES Fermi surface map and constant binding energy contours of the band structure 
of NbAs along the (001) direction at various energies. The square BZ and C 4 rotational symmetry 
at higher binding energies both clearly indicate that the sample is cleaved on the ( 001 ) surface. 
The C 4 violation by certain bands at shallow binding energies is consistent with C 4 screw axis 
symmetry broken by the ( 001 ) surface and clearly shows that the C 4 asymmetric states are 
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FIG. 3: surface states, b, High-resolution ARPES Fermi surface map, with the soft X-ray ARPES 
map from 2c overlaid on top of it to scale (grayscale region), showing the relative positions of 
the Fermi arcs and the Weyl nodes, c, Schematics of the locations of the W2 Weyl nodes in the 
bulk BZ of NbAs and the unusual topology of the Fermi arcs on the (001) surface Fermi surface 
in this compound. The pseudo-spin texture near the Weyl nodes with positive and negative chiral 
charges resembles the magnetic field around magnetic monopoles and antimonopoles, d, ARPES 
Fermi surface map of the tadpole Fermi arcs along P — A and the corresponding Weyl points, e, 
First-principles calculated Fermi surface map of the tadpole Fermi arcs along P — A which shows 
the qualitative agreement with the ARPES data, f, Schematics of the connectivity of the Fermi 
arcs in the first BZ. g, ARPES Fermi surface, and a constant binding energy contour at binding 
energy of 120 meV, of the Fermi arcs along P — A. 
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FIG. 4: Visualization of the co-propagating chiral modes and the structure of the 
tadpole Fermi arcs in NbAs. a, Bottom panel: a circular pipe through the bulk BZ of NbAs, 
enclosing two Weyl points of chiral charge +1. The pipe extends through the entire bulk BZ in the 
kz direction, forming a closed manifold. Because it encloses two Weyl points of chiral charge +1, 
the Chern number on this manifold is +2. Middle panel: by cleaving an NbAs sample on the (001) 
surface, we introduce a boundary on this manifold corresponding to a circle in the surface BZ. Top 
panel: the Chern number of +2 requires that the band structure on this circle in the surface BZ 
have two gapless co-propagating chiral edge modes, b, Our direct experimental observation by 
ARPES of the way in which Fermi arcs connect Weyl points on the (001) surface of NbAs places 
strong constraints on the dispersion of the Fermi arcs. To be consistent with a chiral charge of 
±2 for the projections of the Weyl points on the surface BZ, the tadpole Fermi arcs must disperse 
as shown by the red and blue sheets. The Weyl cones are shown in green, c, Constant-energy 
contours of the tadpole Fermi arc surface states at binding energies e above the energy of the 
Weyl points ew, near the energy of the Weyl points and below the energy of the Weyl points. For 
the Fermi arcs to be co-propagating, the red sheet must grow with binding energy, while the blue 
sheets must disperse toward each other. 














